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A system designed for the calculation of control parameters for a computer based inventory management system according to the 
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System and method for calculation of controlling parameters for a computer based 
inventory management system. 

5 

Technical Field 

This invention relates to inventory management systems and the parameters used therein. In 
10 particular the invention deals with the calculation of optimum reorder parameters, optimizing 
service levels as well as minimizing cost. 

Inventory management systems using suitably programmed computers are not 
unknown as such. As an example may be cited the program system TRITON as supplied by 
15 Baan BV, The Netherlands. This program runs on a computer consisting of at least one 
processor, memory, input and output. 

A theoretical exposition of the functionality of such a system can be found in R.H. Ballou, 
Business Logistics Management, Prentice Hall Inc., 1992 or E.A. Silver and R. Peterson, 
Decision Systems for Inventory Management and Production Planning, John Wiley & Sons, 
20 Inc., 1985 

Background to the invention 

The main aim of inventory management systems is to administrate and maintain a stock of 
25 items, such that ideally any order for any item can be filled from stock, while at the same time 

the stock of items is kept as low as is possible. Constraints such as maximum investment- and 

minimum service levels must be maintained. Excess of stock should be avoided. In other 

words, good inventory management involves providing a high product availability or item 

service level at a reasonable cost. 
30 The service level is defined as the ratio between the number of demands directly served from 

stock and the total number of demands or, alternatively, as ratio of served quantity and the 

total quantity demanded. 
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SL = Served quantity / Total quantity 

The service level is a direct consequence of the, in time on the right moment, reordering 
process for the items in stock. For each item, dependent on the particular supplier of the item, 
5 an estimate of the lead time L may be derived, 
hi order to enable demands to be served from stock during the lead time L, items have to be 
reordered when a sufficient stock quantity, but no more than that, is still available. This way 
the demand during the lead time L can be satisfied. 

To this end items are reordered when the stock quantity reaches a certain level, the reorder 
10 point ROP (see figure 1). This reorder point is chosen to enable demands during the lead time 
to be served, with a certain pre-defined desu-ed service level SLcon. 

Another way of managing the stock consists of inspecting the quantity in stock at fixed time 
intervals (see figure 2.). Reordering occurs to a certain predefined, calculated maximum stock 
level MSL aimed at covering the demand for items during the time intervals between the 
15 inspections, with a certain pre-defined desired service level SLcon. 

Both methods may be used as well, if the items are not available from a supplier, but are 
produced in a batch-wise production with a limited production time and production setup time. 

Management methods are known, however they differ from the invention in their purpose, 
20 their implementation and their mathematical basis. 

Patent 5,287,267, Feb. 15,1994 

This patent describes a method of controlling an mventory of parts, in particular components 
needed for the fabrication of products, the same part sometunes being used in a plurality of 
25 products, where the actual demand for the products is unknown. Based on estimated forecasts 
and confidence intervals for the different products, and bills of material, vvhich express the 
quantities of the parts which are needed as components for each of the products, estimates of 
the total number of parts required are obtained. 

The optunum to an objective aimed at minimization of excess stock, is found using a standard 
30 iterative procedure where the solution is found by use of a parametric search on the value of a 
Lagrange multiplier. 
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This patent describes a method of managing an inventory of fast moving consumables, 
predicting near future consumption by the use of statistics based on the recent past, sampling 
daily consumption patterns, differentiated by day-types. The method is heuristic and uses input 
and control by human expert knowledge. 

5 

Patent EP 0 733 986 A2, Sept. 25, 1996 

This patent describes a method of optimizing an inventory, based on a number of different 
criteria, such as a selected inventory investment or a service level constraint. A number of 
initial parameters such as forecasts are predefined. A level of an inventory investment or an 
10 inventory level constraint is found by marginal analysis. Gradients or slopes of constraint 
functions are used to improve the fit of an ensemble of parameters for the items concemed, to 
predefined service levels or investment constraints. Step-wise improved (on average) values 
for the ensemble can be obtained. This method is a refinement method for an initially known 
ensemble of parameters where optimal individual results are not guaranteed. 

15 

Aim of the invention 

Current methods of calculating the reorder point and maximum stock level MSL can be 
20 improved, in particular by improving the accuracy of the estimates of the variance of the 
consumption over a certain period of time. A mathematical method avoiding sampling by 
periods of time, deriving the solution by solving a mathematical expression, using the 
historical data directly, yields the information required. The invention indicates how this aim 
can be realized. 

25 

Short description of the invention 

In agreement with the aim given the present invention provides a system and method to 
calculate optimal inventory control parameters for the use in a computer based inventory 
30 management system, which comprises at least one processor, memory, input and output. The 
system consists of a computer program, which cMculates the moment of feoMering^ a^ 
the quantity to reorder. The calculation is characterized by the use of moments of the single 
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demand probability functions P(Q) for each item in the inventory to cdculate the reorder point 
ROP and reorder quantity or tiie maximum stock level MSL using the formulae 

E{R) = AlxE(L)xE(Q) 

5 

And 

Var(R) = A2 x E(L) x Var(Q) + A3 x E(L) x E(Q)2 + AA^ x E(Q)2 x Var(L) 
10 Where 



Aj parameter depending on the distribution function of Tj 

Tj Time elapsed between the demands j-1 and j. 

L The lead time, reorder period or production time 

15 R Consumption in units over L 

Q(i) Quantity in the single demand i during lead time 

E(R) Expected consumption of item over L 

E(L) Expected value of L 

E(Q) Expected value for Q for any request i 

20 Var(R) Variance of R 

Var(Q) Variance of Q 

Var(L) Variance of L 



An alternative embodiment of the system comprises using the following formulae 

25 

n 

P(R)= S wjxPj(R) 
j=0 

30 Where 

P(R) the probability distribution of R 

P(Q) the probability distribution of Q for any request i 
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wj 



Pj(R) 



the jth selfconvolution of P(Q) (the joint probability for the total quantity of j 
requests) 

the statistical weight of the corresponding joint probability distribution for 1 .... 
n simultaneous demands in lead time L. 



5 



The reorder quantity is calculated from: 



ROP 



SL(ROP) = |P(R)dR 



10 



0 



In yet another alternative embodiment the following formula is applied 



15 SL(ROP)= 2 wjxFj(ROP) 
j=l 

Where 

20 SL(ROP) Service level as function of the ROP 
Fj(ROP) The cumulative distribution function of Rj 

In still another alternative embodiment a completely general formula is applied: 

25 n 

L(R) = 2 wjx[L(Q)]j 
j=l 



n 



Where 



30 



L(Q) 
L(R) 



the Fourier transform of P(Q) 
the Fourier transform of P(R) 
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P(R) is found by back-transforming the result just once. 
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P(R) = Fourier transform of ( L(R) ) 

5 From this expression ROP, given a desired service level SLcon, can be obtained directly. 

Brief description of the drawings 

Figure 1 depicts the Q-system. If the stock reaches the reorder point ROP, an order is 
10 created. The reorder point ROP is calculated in such a way as to enable serving demands for 
items during the lead time L, with a certain pre-defined desired service level SLcon. 

Figure 2 depicts the fixed order cycle method, or the P-system, where on periodic time- 
intervals Tp an order is created. At periodic intervals, the quantity in stock is reviewed, and an 
order is created to re-supply the stock. The maximum stock level MSL enables demands to be 
15 served between inspections to a certain predefined service level SLcon. The lead time L is 
taken into account as well. 

Figure 3 depicts the relation between the stock, the lead time and variability of the lead 
time, and demand and variability in demand. The in reality step-wise decrease of the inventory 
stock is depicted as straight lines, the slope representing the total quantity demanded in the 
20 lead time. 

Figure 4 depicts a probability distribution which can be used to determine the reorder 

point. 

Figure 5 depicts two ways of representing the historical consumption over time of an 
item: First, each demand for an item is represented by a bar, the length of which is 
25 proportional to the quantity requested. Alternatively, the historical consumption is represented 
by sampling the total quantity demanded over a specific period of time. 

Figure 6 depicts two different demand pattems A and B, each are represented in two 
ways according to figure 5, as single demanded quantities representations Al and Bl, and as 
the demanded quantities sampled over periods, A2 and B2. 
30 Figure 7 depicts the same demand pattern A, in two ways A2 and A3: As sampled by 

periods using tiie same sampling period intervals, the sampling period intervals of A3 have 
been shifted with respect to A2. 

Figure 8 depicts a number of graphs clarifying the detailed description of the invention. 
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Figure 9 depicts a diagram illustrating an aspect in the detailed description of the 
invention. 

5 Detailed description of the invention 

Using the values of the parameters calculated by the control system according to the 
invention, allows the inventory management system to optimize the service levels while 
minimizing the stock inventory and avoiding excess stock. The logistical parameters 
10 calculated by the system are: 

For the Q-system (R.H. Ballou, Business Logistics Management, Prentice Hall Inc., 1992) the 
reorder point and the reorder quantity. 

For the P-system (R.H. Ballou) or fixed order cycle method, the order to level or maximum 
stock level MSL and the periodic review period Tp. 
15 For both Q-system and P-system the break quantity or exceptional demanded quantity 
threshold and lead times are calculated as well. 

Currently, the reorder point is calculated as expected demand E(R), the so-called stock reserve 
SR, to which a quantity is added, the so-called safety stock SS, to cover larger than expected 
demand, caused by the variability in demand and lead times. See figure 3. 

20 The reorder point ROP for Q-systems and order to level MSL for P-systems are calculated in 
exactly the same manner, the only difference being that ROP depends on the lead time L, 
whereas MSL depends on Tp + L. Therefore the equations below are given for ROP only. 
Similar calculations are equally applicable if the items are not re-supplied by ordering from a 
supplier, but are produced in batch-wise production with a limited production time and 

25 production setup time. 

Current methods of determining reorder points generally use the following mechanism to 
estimate ROP: The probability distribution of P(R) is assumed to be a normal distribution 
(Figure 4). If demand R > ROP, the demand will not be served. This constitutes the area or 
fraction p. The estimate of ROP is given by 

30 

ROP = E(R) + ZxSd(R) 
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Where Z is given by Z = G-l(l-p) resulting in a service level s = (1-p) x 100% (G is the 
standard normal distribution function). 

Sd(R) is estimated by several methods, the one most commonly used is based on the mean 
weighted average absolute difference of the consumption and forecast over a number of 
5 historical periods: 



n 

MAD = S wi X I Forecast (period i) - Real demand (period i) | 
10 i=l 

Where the weights wi correspond generally to exponential smoothing, the factor 1.25 is used 
to convert the mean-absolute-deviation-value into the theoretical desired square-root-mean- 
squared-deviation 

15 

Sd(R) = 1.25 X MAD 

The value of MAD should be corrected for the difference between L and the forecast time on 
which MAD is based. 

20 A second method is based on an estimation of Sd(R) by direct estimation of the variance of the 
demand sampled in periodic intervals, to which an additional term is added, allowing for the 
variance of the lead time. 

Var(R) = E(L) X Var(D) + [E(D)] ^ X Var(L) 

25 

Both of these methods are based on sampling the total quantity demanded D by periodic 
intervals only, rather than using the statistics based on individual demands, hi this context 
terms such as "variability in demand" and "variance of demands" are not uncommon, however, 
this does not imply the use of individual demand statistics. 
30 Note that sampling of demands in intervals actually leads to loss of information: 

Li figure 5, a bar represents a demand for a quantity q, the length of the bar is proportional to 
the quantity. A box represents the total demand D within a given period. 
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Two different patterns of demand are represented in figure 6 in pattern Al and Bl 
respectively, together with the presentation of the same demand patterns sampled in periodic 
intervals in A2 and B2. 

5 Clearly, sampling in a periodic way leads to the false conclusion that the non-identical patterns 
shown are in fact identical. 

Moreover, as is depicted in figure 7, small changes in - or shifts of the period of sampling of 
the same pattern may lead to radically different estimates for the variance of the distribution as 
10 can be seen from A2 and A3, both periodic presentations of the same demand pattern Al . 
(Note that where the variance is highly dependent on the way the pattern is sampled, changing 
the sampling or shifting the presentation does not influence the mean or average expected 
quantity for a sampling period.) 

15 Only the service level during lead time is important, the service level outside lead time is 
always at least equal to this. In the context of the invention, ROP is estimated from historic 
data while avoiding the problems mentioned above. 
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P(R) is derived from P(Q), where P(Q) is the empirical probability distribution, based on 
historical frequency data (figure 8.) which may be time-weighted using an empirical weighting 
5 scheme wq, down weighting inaccurate data and data obtained from a long time in the past. 
If the weighted frequency distribution is used, the sum of the corresponding weights is taken 
for all demanded quantities pertaining to a certain interval, mstead of the number of times a 
demanded quantity pertains to the same interval. 

10 The historical frequency distribution is obtained by sampling Qi on a number of intervals ql, 
q2 ... qm and counting the frequency of occurrence for each interval. 

F(qm) = frequency Qi in interval qm = sum of 1 (Qi in interval qm) 

1 5 The weighted frequency distribution is obtained by adding the to Qi corresponding weight wq 
instead of adding 1, when summing the number of occurrences for each interval. 

Fw(qm) ~ sum of wq (Qi in interval qm) 

20 

Given a number of demands, for instance 3 (figure 9.), the joint probability distribution P3(R) 
for the total consumption R of 3 demands, can be constructed from the single demand 
probability distribution P(Q). 

25 The demand in the lead time, R, may arise from different numbers of demands j, each with a 
certain probability of occurring wj, and a probability distribution Pj(R) for the total quantity of 
thej demands. 

Therefore P(R) is derived as a series. 



30 n 

P(R)= S wjxPj(R) 
(formula 1) 
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Where Pj(R) is the jth selfconvolution of P(Q) and wj is the statistical weight of the 
corresponding joint probability function for 1 ... n simultaneous demands in the lead time L. 
In practice values of n>100 need not be considered, as then the alternative approach given 
5 below is valid. 

Under the assumption of the distribution of Ti being known - e.g. an exponential or truncated 
normal or Weibull, etc. - the coefficients wj can be calculated directly. This calculation of the 
coefficients can be modified to include the function P(L). 
10 In the case of the distribution of Ti being exponential, for wj a Poisson distribution is the 
result, with a mean demand number density A equal to the expected number of demands 
during L divided by L. Calculation of P(R) is then straightforward if the functions Pj are 
known. 



15 The functions Pj are not easily obtained directly, however estimates of sufficient accuracy of 
these convolutions are obtained by Fourier transformation of P(Q) and back-transforming the 
jth power of the transform obtained. However, P(R) may be obtained directly by first summing 
over j the jth powers of the transform of P(Q), 

20 n 

L(R) - S wj X [L(Q)] i (formula 2) 

using the weights wj, and back-transforming the result just once, 

25 

P(R) = Fourier transform ( L(R) ) (formula 3) 



From this expression, ROP, given a desired service level SLcon, can be calculated easily. The 
reorder quantity is now calculated in a conventional manner. 

30 

ROP 

SL(ROP) = I P(R) dR (formula 4) 

0 
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In order to improve the accuracy and avoid series termination effects in Fourier space, the 
Fourier-transform L(Q) can be multiplied with the, in the reciprocal space defined, weighting 
factor wL. 

5 

n • 
L(R) = E wj X [ wL X L(Q)] J (formula 
2B) 

10 

Alternatively, for certain probability distributions of P(Q), such as a normal or a gamma 
distribution, Pj(R) can be found analytically for all values of j. The ROP can be calculated 
using formula 4 directly or from the cumulative distribution functions Fj of Pj(R) in formula 5, 

15 

n 

SL(ROP) = E wj X Fj(ROP) (formula 5) 

j=l 

20 If the number of demands during lead time is sufficiently large, P(R) can be considered to be a 
normal distribution, hi this case P(R) is fully defmed by its mean value and variance. The 
mean may be obtained from 

E(R) =AxE(L)xE(Q) (formula 6) 

25 

And the variance may be obtained from 

Var(R) = A x E(L) x Var(Q) + A x E(L) x E(Q)2 + A^ x E(Q)2 x Var(L) (formula 7) 

30 The advantage is obvious in terms of speed and efficiency. However, it must be stressed that 
this approach is only valid if de number of demands is high. 

Note that formula 7 comprises the (first and second moment, mean and variance) moments of 
the single demanded quantity probability distribution. The E(Q), Var(Q) may be time- 
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weighted using an empirical weighting scheme wq, down weighting inaccurate data and data 
obtained from sampling a long time in the past. 

E(Q) = E wqxQj/2wq 

5 

And 

Var(Q) = 2 wq X (Qj - E{Q))^ I S wq 

10 If in the Q-system the quantity in stock drops below the reorder point ROP by serving a 
demand, the calculation does not start at the reorder point ROP but at a point well below the 
reorder point. In order to compensate for this effect a correction term may be applied to the 
above mentioned formulae. 

1 5 Variables used in the description, for each item in the inventory: 

MAD Mean average absolute deviation between forecasted and actual consumption over a 
number of historical periods 

D Demand rate, total quantity demanded by period 
20 Var(D) Variance of D 

L The lead time, re-order period or production time 

P(L) Probability distribution function of L 

E(L) Expected value of L 

Var(L) Variance of L 
25 R Consumption in units over L 

P(R) Probability distribution function of R 

E(R) Expected consumption of item over L 

Var(R) Variance of R 

Sd(R) Standard deviation of R 
30 Tj Time elapsed between the Q-l )-th and the j-th demand 

P(T) Probability distribution function of Tj 

Rj Consumption in units of j demands over L 

Pj(R) Probability distribution function of Rj 
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wj the statistical weight of the corresponding joint probability function Pj(R) for L.j 
simultaneous demands in lead time L. 

Fj (ROP) The cumulative distribution function of Rj 

Q(i) Quantity in the i-th single demand during lead time 
5 wq Empirical weighting scheme, down weighting inaccurate data 

P(Q) Probability distribution function of Q for any i 

E{Q) Expected value for Q for any i 

Var(Q) Variance of Q 

A Mean demand number density 
10 L(Q) Fourier transform of P(Q) 

L(R) Fourier transform of P(R) 

wL Weighting factor in reciprocal space 

ROP Reorder point 

SL Service level 
15 SL(ROP) Service level as function of the ROP 

SLcon Desired service level 

Tp The time between periodic reviews 

MSL Order-to-level or maximum-stock-level 



20 
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Claims: 



1 . System designed for the calculation of optimal control parameters for use in a computer 
based inventory management system comprising a computer program designed to calculate the 
5 reorder point or the maximum stock level MSL, characterized by the use of the single 
demanded quantity distributions for each item in the inventory in order to determine the 
reorder point or, depending on the reorder policy, the maximum stock level MSL for periodic 
reviews. 

10 2. System designed for the calculation of optimal control parameters for use in a computer 
based inventory management system, according to claim 1, characterized by the use of the 
smgle demanded quantity distributions for each item in the mventory in order to determine the 
reorder point or depending on the re-order policy, the maximum stock level MSL for periodic 
reviews, using the formulae 



15 



E(R) =AlxE(L)xE(Q) 



and 



20 Var(R) = A2 x E(L) x Var(Q) + A3 x E(L) x E(Q) ^ + A4 ^ x E(Q) ^ x Var(L) 



where 



Aj 
25 Tj 



30 E(L) 



L 



R 



Q(i) 
E(R) 



E(Q) 

Var(R) 

Var(Q) 



parameter depending on the distribution function of Tj 

Time elapsed between the (j-l)-th and the j-th demand 

The lead time, re-order period or production time 

Consumption in units over L 

Quantity in the i-th single demand during lead time 

Expected consumption of item over L 

Expected value of L 

Expected value for Q for any request i 



Variance of R 



Variance of Q 
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Var(L) 



Variance ofL 



3. System according to claim 2 characterized in that the distribution Tj is exponential in 
which case Al = A2 = A3 = A4 

5 

4. System according to claim 2 or 3 characterized in that it uses weights wq to improve the 
accuracy according 

E(Q) = SwqxQj/2wq 

10 



Var(Q) = E wq X (Qj - E(Q))2 / Z wq 

15 

5. System designed for the calculation of optimal control parameters for use in a computer 
based inventory management system, according to claim 1, characterized by the use of the 
single demanded quantity distributions for each item in the inventory in order to determine the 
reorder point or depending on the reorder policy, the maximum stock level MSL for periodic 
20 reviews, using the following formulae 



and 



n 



P(R)= E wjxPj(R) 
j=0 



25 



where 



P(R) 
P(Q) 
30 Pj(R) 



the probability function of R 

the probability function of Q for any request i 

the jth selfconvolution of P(Q) (the joint probability for the total quantity of j 
requests) 

the statistical weight of the conesponding joint probability function for 1 .... n 
simultaneous demands in lead time L. 
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after which the reorder quantity is now calculated as: 
ROP 

5 SL(ROP) = J P(R)dR 
0 

6. System designed for the calculation of optimal control parameters for use in a computer 
based inventory management system, according to claim 1, characterized by the use of the 

10 single demanded quantity distributions for each item in the inventory in order to determine the 
reorder point or depending on the reorder policy, the maximum stock level MSL for periodic 
reviews, using the formulae 

n 

15 SL(ROP)= S wjxFj(ROP) 

ri 

where 

20 SL(ROP) Service level as function of the ROP 
Fj(ROP) The cumulative distribution fimction of Rj 

7. System designed for the calculation of optimal control parameters for use in a computer 
based inventory management system, according to claim 1, characterized by the use of the 

25 single demanded quantity distribution for each item in the mventory in order to determine the 
reorder point or depending on the reorder policy, the maximum stock level MSL for periodic 
reviews, using the formulae 

n 

30 L(R) = I wjx[L(Q)]j 
j=l 

where 
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L(Q) the Fourier transform of P(Q) 

L(R) the Fourier transform of P(R) 

5 after which P(R) is found back-transforming the result just once: 

P(R) = Fourier transform ( L(R) ) 

whereafter the reorder quantity is calculated as: 

10 

ROP 

SL(ROP) = J P(R) dR 
0 

1 5 8. System designed for the calculation of optimal control parameters for use in a computer 
based inventory management system, according to claim 7, characterized by the use of a 
weighting factor wL: 

n 

20 L(R) = E wjx[wLxL(Q)]-i 
j=l 

where 

25 L(Q) the Fourier transform of P(Q) 

L(R) the Fourier transform of P(R) 

wL weighting-factor in the reciprocal space 

after which P(R) is found back-transforming the result just once: 

30 

P(R) = Fourier transform ( L(R) ) 
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Figure 3 
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Figure 4 
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Figure 5 



A demand for a quantity q is represented by a bar, the length of the bar is proportional 
to the quantity. The total demand D within a given period is represented by a box. 



Figure 6 

Two different patterns of demand Al and Bl are represented below, together with the 
presentation of the demand sampled in periodic intervals in A2 and B2. 
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Figure 7 

A small shift of the period of sampling of the same pattern may lead to radically 
di£ferent estimates 

A2 
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